AN APPLICATION OF THE MAX-PLUS SPECTRAL THEORY 
TO AN ULTRADISCRETE ANALOGUE OF THE LAX PAIR 
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Abstract. We study the ultradiscrete analogue of Lax pair proposed by 
Willox et al.|lll. This "pair" is a max-plus linear system comprising four 
equations. Our starting point is to treat this system as a combination of two 
max-plus eigenproblems, with two additional constraints. Though infinite- 
dimensional, these two eigenproblems can be treated by means of the "stan- 
dard" max-plus spectral theory. In particular, any solution to the system can 
be described as a max-linear combination of fundamental eigenvectors associ- 
ated with each soliton. We then describe the operation of undressing using 
pairs of fundamental eigenvectors. We also study the solvability of the com- 
plete system of four equations as proposed by Willox et al. 11. . 



1. Introduction 

1.1. Motivations and purposes. We consider the system of four equations 
max($[;\ - k, = + max(t/(!\ - 1, 

max($[;V^ - k, = ^i^'^ + niax(C//*) - 1, -[//!\), 



(1.1 



niax($p^ +k-uj, + + fc - 1) = $[*+^^ 

max(<r\ <\+c//*'-i)=<i>r^ 



whiclr appeared in tlie work of Willox et al. [TT] . Here we assume that the potential 
is known, and that and the solutions satisfy the conditions 

(AU) and (A$) written below, see Subsection 12.21 

System p.ip plays the role of the Lax pair for the ultradiscrete KdV equation 



(1.2) [//*+^) = min(l - C//*), ^ Ui'^ ~ ui'^\ 

k— — oo 

which describes the dynamics of Box& Ball system of Takahashi and Satsuma [7] . 
Willox et al. [11] show how solving (|l.ip helps to calculate the phase-shifts of 
solitons after interaction in the case of the real initial U and, more generally, to 
solve equation ()1.2|) at all times. 

Very briefly, the relation of ()1.2|) to the classical discrete and continuous KdV 
equations is as follows. It was shown by Tokihiro et al. liSl that equation (|1.2|) can 
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be obtained as ultradiscrete limit (or Maslov dequantization) of the discrete KdV 
equation 

(1.3) ^_J_^j(,W s>0 

written by Tsujimoto and Hirota [9 . This equation turns into the famous Lotka- 
Voherra equation by taking the continuous hmit (5 0) [B] , and the Lotka- Voherra 
equation is also known as an integrable discretization of the classical KdV equation. 
See |H] and [TT] for more explanation. 

The intention of this paper is to build a max-plus linear theory of (|l.ip . To our 
point of view, such theory is lacking in [TT] , where it is claimed that system (II. ip 
is always solvable, but without going into the details of the proof. As we will see, 
the theory of system (jl.ip is nontrivial and to the author's knowledge this kind of 
problems never appeared in the max-plus literature and could be of its own interest. 
Namely, we have two infinite max-plus eigenproblems represented by the first two 
equations of (where the eigenvalue is necessarily 0), and two connections 

between them represented by the last two equations. Thus we are led to study two 
related (but different) eigenproblems at the same time, taking into account some 
additional constraints. 

In this paper we do not address the solvability statement of [11 in full strength, 
and rather concentrate on developing the spectral theory associated with the first 
two equations of (jl.ip . These are two closely related infinite max-plus eigenprob- 
lems of a special kind. The theory of such problems was developed by Akian, 
Gaubert and Walsh pi, and it could be applied here. However, we notice that as- 
suming conditions (AC/) and (A$) on U and $, the problem can be reduced to the 
more usual finite max-plus spectral theory as described in the monographs [3l SI [6] . 
Namely with each soliton of C/'^*^ we can associate a pair of fundamental eigen- 
vectors, and any solution of the first and of the second equation of (11.11) appears 
as their max-plus linear combination, see Proposition 12.71 and Theorem 12.81 Thus 
we describe the set of all solutions to the first two equations of (|l.ip with natural 
asymptotic behaviour. 

Next we consider the procedure of undressing the initial potential U by means of 
a pair of fundamental eigenvectors. It follows that in a natural special case when the 
"interior" of a soliton in U consists of I's, this soliton disappears after undressing, 
and the rest of the potential gets shifted by one position towards the soliton. Note 
that in the undressing procedure of [TT], it is demonstrated that the other solitons 
may change their form. This effect does not happen in our case, since we use the 
fundamental eigenvectors. 

Finally we treat the complete system (|l.ll) . In the case when U has no massive 
solitons (J/f*'* -I- t/f+i < 1 for all i) or when U has just one massive soliton, we 
confirm that (|1.1[) is solvable by showing that any pair of fundamental eigenvectors 
is a solution. In the case of several massive solitons we show that to the contrary, 
no pair of fundamental eigenvectors is a solution, so that a combination of these 
fundamental eigenvectors satisfying (jl.ip has to be guessed. 

1.2. MfLx-plus spectral theory. Algebra max-plus is developed over the real 
numbers R completed by the least element — oo, with arithmetical operations a (g) 
b := a + b ("multiplication") and a® 6 := max(a, b) ("addition"). The new "zero" is 
— oo and the new "unity" is 0. This arithmetics is extended to matrices and vectors 
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in the usual way so that 

{A B),fe = Aij (g) Bjk , {A © B)ij = A^J © B^J , 
j 

for matrices A and B of appropriate sizes. We wiU be interested only in the max- 
plus spectral problem 

(1.4) ^©<I> = A(g)$, 

that is, trying to find for a matrix A e (MU{— oo})"^" a parameter A G MUj— cxo} 
such that there exists a vector $ satisfying (|1.4I) with not all components equal to 
— oo. 

Max-plus spectral theory uses the following graph-theoretical concepts: 

1. Associated graph G{A) = {N, E) with set of nodes N — {1, . . . , n} and set of 
edges E = : Aij ^ — cxd} weighted by w{i,j) = Aij. The concept of weight is 
extended to paths P = {io ^ ii ^ . . . ^ ik) , defining the weight of P by 

w{P) :— Aigi^ ■ . . . ■ Ai^_-^i^ 

Closed paths P having iQ = are called cycles. 

2. Critical graph G'^{A) comprising all nodes and edges that belong to the cycles 
(ii, . . . , ik), on which the maximum in 

/-I r\ \ / 4\ -^1112 © ■ • • © Ai.i, 

(1.5) X{A) = max max — — ; — 

is attained. Such cycles are called critical, and so are all nodes and edges of the 
critical graph. Being made from cycles, the critical graph is completely reducible, 
i.e., it consists of several isolated strongly connected components. 

3. Saturation graph Sat($) consisting of all nodes and edges satisfying Uij+^j = 
A + $i, that is, attaining maximum on the l.h.s. of (jl.4p . 

The following theorem explains some properties of the saturation graph and its 
relation to the critical graph. It is well-known but we give a short proof for the 
reader's convenience. 

Theorem 1.1. Let A e {RU {-oo})"^". Suppose that $ satisfies A © $ = $ and 
has all components finite (i.e., not —oo). Then 

1. Each node has an outgoing edge in Sat($), 

2. Each cycle in Q has total weight not exceeding 0, 

3. The cycles o/Sat($) are precisely the cycles of Q'^{A) 

Proof. 1.: If $ is an eigenvector then for each i there exists j such that Aij+^j = ^i. 
2. and 3.: Let (zi, . . . , jfc) be a cycle in G{A). Then we have 

Combining these inequalities and cancelling $ we get Ai^i.^ © ... © Ai^i^ < 0, which 
shows 2. Note that Ai^i^ + . . .+Ai^i^ < is equivalent to having < ^ii 

for some I, which implies 3. □ 

Note that Theorem 11.11 generalizes to the case when the matrix A is infinite- 
dimensional but each row has a finite number of real entries. This is the case that 
we will have to work with when analyzing (jl.ip . 

For A G (K U {— oo})"^", a formal analogue of (/ — A)^^ can be defined as 

(1.6) A* ^I®A®A^®..., 
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where / is the m£LX-plus identity matrix, having on the diagonal and — oo 
otherwise. A* converges and can be truncated to I Q) A Q) ■ ■ . if and only 

if X{A) < 0. Note that entries of A*, denoted by A*j, have a principal path 
interpretation: for i j this is the greatest weight of a path connecting i to j. 

When X{A) = it can be shown that any column A*^ of A*, whose index i 
belongs to the critical graph (i.e., critical column of A*), is an eigenvector of 
A. Such columns are called the fundamental eigenvectors. The eigenspace of A 
can be described more precisely as follows. 

Theorem 1.2. Let A e (M U {-oo})"^". // X{A) anrf $ satisfies A (g) $ = 
then there exist aiGMUj— oo} such that 

where S G is any index set containing precisely one index from each 

strongly connected component ofG'^iA). 

That is, each eigenvector of ^ is a max-linear combination of the fundamental 
eigenvectors. 

Theorem 11.21 can be found in several monographs on max-plus algebra O |4j [6] . 
The max-plus spectral theory (both finite- and infinite-dimensional) has applica- 
tions ranging from railway scheduling [6] to Frenkel-Kontorova model in solid state 
physics [5], [To], and the crop rotation problem in the agriculture 0. 

2. Simplified Lax pair 

2.1. Solitons and critical graphs. For most of this paper we will consider a 
simplified version of the Lax system for udKdV consisting of the first two equations 
of (|l.ll) . which we rewrite as 

(2.1) max[$(;V^ +6.~k, + S^] - $f 

where7, = min(C/f \ 1 - ) , 6,^ min(C/f!\ A~U^'^)- 

Further we will fix t and denote u, := f/f \ $(1) := and $(2) $(*+!) ^ 

We distinguish between two cases: 
(CI) When v^up = supj(ui + u^+i) < 1. 
(C2) When v^up = supj(ui -I- u^+i) > 1. 

Note that the borderline case supj Ui + u^+i = 1 can be regarded in both ways, 
not leading to any contradiction. 

We observe that (|2.ip is a combination of two max-plus eigenproblems ^(7) (S) 
_ (f^].g^ equation) and A{6) (g) = <i>(^) (second equation), where the 
coefficients of ^4(7) and A{S) can be written as follows. 

Lemma 2.1. In the case (CI), 

^2 ^(7)1+1, i = "i+i, ^(7)^^1+1 = Ui - fc, 

A{6)i+i,i = Uj, y4((5)i,i+i = - k. 
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S{S): 



U: 




Ul Ui + i Ul+2 

U;-2 — k lt;_i — fc Ul — k 




Figure 1. Case (CI): a fragment of the associated digraphs 



Lemma 2.2. In the case of (C2) 

^(7)^+l,^ 



Ui+i, if Ui + Ui+i < 1, 

1 - Ui, if Ui + Ui+i > 1 



(2.3) 



Ui - k, 

1 - Ui-i 



if Ui 
k, if Ui 



Ui-l < 1, 
Ul-l > 1 



u 

1 — u 



if Ui 
if Ui 



Ui+l < I, 
Ui+l > 1 



Ui_i — k, if Ui + Ui_i < 1, 
1 — Ui — k, if Ui + Ui_i > 1 



The proofs are straightforward. We proceed with the following crucial definition. 

Definition 2.3 (Solitons). 1. In the case of (CI), soliton is a sequence of 
indices {1,1 + 1, . . . ,1 + s) such that ui + = uj+s-i + = Wsup, 

while ui-i +ui < Vsup and ui+s + ui+s+i < Wsup- 
2. In the case of (C2), soliton is a sequence of indices {1,1 + 1, . . . ,1 + s) 
such that Ul + u;_|_i > 1, . . . , ii;+s_i + Ui_|_s > 1, while ui^i + u; < 1 and 

Ul+s + Ui+s+l < 1- 



If in the equation of cellular automaton (|1.2I) we assume that Ui*^ = Ul'''^^' = 
for all i < —N where N is sufhciently large (see condition {AU) below), then 
its dynamics can be computed explicitly. In the case (CI) the whole vector t/*^*^ 
gets shifted by one position to the right. In the case (C2) the behaviour is more 
complex. Like in the classical theory of KdV, the solitons (as defined above) move 
with different speed depending on their mass (not defined here). After interaction 
they emerge again with a phase-shift, as described by Willox et al. [TT]. So it can be 
argued that Definition 12.31 has a "physical sense" only in the case (C2). However, 
as we show below, the theory of eigenproblems (12.11) is similar in both cases. 

Consider, with Lemmas 12.11 and 12.21 in mind, the associated weighted digraphs 
^(7) and G{S) of matrices ^(7) and A{6). We are going to study the critical cycles, 
i.e., the two-cycles with the greatest total weight, and the critical graph, consisting 
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Ui^l — k Ul — k 1— U/ — fc 




Ul 1 - Ul 1 - Ui + i 

M;_2 - k ui^i - k 1 - ui+i - k 



U: 




Ul-l 1 - Ul + i 1 - Ui+2 

Ul-1 Ul Ui + i Ui+2 



Figure 2. Case (C2): tail of a soliton. Cycle (2,3) of Q{j) is critical. 

1 - Ul+s-3 - k 1 ~ Ui+s-2 - fc 1 - ui+s-i - k 
Q{,l) : 1 2 3 4 




1 - U;+s-2 1 - Ui+s+l 

1 - Ui+s-2 - fc 1 - - fc 1 - Ui+s - k 




Figure 3. Case (C2): head of a soliton. Cycle (3,4) of Q{5) is critical 

of all nodes and edges on the critical cycles. We relate the strongly connected 
components of critical graphs to solitons, and we give a formula for the greatest 
total weight when the solitons exist. 

The case (CI) is displayed on Figure [TJ Clearly, solitons correspond to the 
strongly connected components of the critical graph (if it is non-empty) , consisting 
of the two-cycles with the greatest total weight maxi(ui -I- u^+i — k). 

In the case of (C2), we give only fragments of these digraphs corresponding to the 
tail (i.e., the left end) and the head (i.e., the right end) of any soliton. The reader 
may assume fc = 1, which will follow from ProDOsition l2.61 under some assumptions 
on U and <&. See Figures [2] and [3l 

All simple cycles of these graphs have length two. If we are in the case (C2), 
then the critical cycles of ^(7), i.e., the ones with the greatest sum of the weight 
of edges equal to 1 — fc, are in the tails of all massive solitons, between the nodes 
marked as 2 and 3. Likewise for Q{5) the critical cycles (of the same total weight 
1 — fc) are in the heads of all massive solitons, between the nodes marked as 3 and 
4. Other cycles in the soliton can be also critical, if and only if ui -\- Wi+i = 1 for 
the corresponding i. We obtain that the greatest total weight of two-cycles is 1 — k. 
What we observed so far can be summarized as below. 
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Theorem 2.4. // the solitons exist, then the greatest total weight of two-cycles is 

(2.4) min(max(Mi_i -f Mi), 1) — fc. 

i 

In this case the critical graphs of both A{'j) and A(6) are non-empty, and their 
strongly connected components are in one-to-one correspondence with the solitons. 

2.2. Reduction to the niELx-plus spectral theory. Further we will assume the 
following requirements on the potential U and on the solution $. 

(AJ7) There exists such that Ui = for all z > iV and i < -N. 

(A$) There exist arbitrarily large N' and N" such that $Ar/ = <i>Ar/+i and 

^-N"-l = $-7V" - li- 
lt will be clear that (A^) is equivalent to the bound state condition of [llj : 
that tends to — oo when i ^ —oo, and that <i>i is constant for i > N' for some 
N'. 

In what follows we treat ^(7) and A{5) at the same time, denoting them by A. 
The associated digraph will be denoted by Q. We will need the following immediate 
observation (following, for instance, from Lemmas 12.11 and 12.21) . 

Lemma 2.5. We have Ai.j+i = —k and ^i+i.i = for i > N and i < —N — 1. All 
weights of cycles (i, i + 1) equal — fc. These cycles have weight (all of them) if and 
only ii k — 0. 

With assumption (AU) we observe that Q always has cycles attaining the greatest 
total weight, and that this weight is given by ()2.4|) . Indeed, using Assumption (AU) 
we see that m-i^Ui = for all i ^ [— A^+l, A^] . In the case when there is i such that 
Ui-i + Mi > 0, this guarantees existence of solitons and leads to <\2A\ . Otherwise 
the claim is trivial with the greatest weight equal to — fc (see Lemma 1^75)1 . 

The saturation digraph of $ can be introduced as in the introduction: 

(2.5) (z,j)eSat($)^^,, 

We now proceed with a proof (elementary but tedious) that with assumptions 
(At/) and (A<I>), the solutions $ are essentially the eigenvectors of the submatrix 
extracted from the interval [— A^ — 1, A^ + 1]. In Proposition 12.61 we will show that 
Sat(<I>) has an outgoing edge from all nodes in [— A^ — 1, A^ + 1] pointing inside 
this interval. We will confirm that the edges of Sat($) outside [—N — 1, A' + 1] 
are directed to this interval, and that the formula for k is as suggested by (|2.4p . 
Based on these observations on Sat($), we show in Proposition l2 . 71 that any solution 

can be obtained as a unique extension of an eigenvector of the submatrix of A 
extracted from [— A^ — 1, A^+1]. The description of solutions in terms of fundamental 
eigenvectors is obtained in Theorem 12.81 

Proposition 2.6. Suppose that U satisfies assumption (At/), $ satisfies A®^ — ^ 
and assumption (A$). Then 

1. Sat($) contains all backward edges i — > i — 1 for i > A^ and all forward 
edges i — i + 1 for i < — A^ 

2. In the restriction of Sat(<l') to [— A^ — l,N + l] for / > 1, every node has an 
outgoing edge. 

3. fc is given by 

(2.6) fc = min(max(ui_i + Ui), 1). 
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Proof. As (A<i>) and {AU) are satisfied, there exist arbitrarily large N',N" > N, 
such that 

^,-1 = -/c + $i = Ai^i,i + for i = -N", 

2.7 

= $i = Ai+i^,, + for i ^ N'. 

Thus the edges iV' + 1 ^- TV' and -TV" - 1 ^ -N" belong to Sat($). 

First we have to treat the exceptional case k = 0. In this case the cycles (iV', iV'+ 
1) and {-N", ~N" - 1) belong to Sat($). Hence by Theorem O they have the 
greatest cycle weight in 5, which is equal to 0. In this case also all cycles for 
(i — 1, i) for i < —N and {i, i + ior i > N have weight and belong to Sat(<I>) by 
Theorem ll.li hence part 1. Having Ui + u^-i > is impossible in this case, as it 
leads to cycles with a positive weight. If Ui + Ui-i < for all i, then the greatest 
cycle mean is max(Mi + = 0, which equals 0. Hence also part 3. For part 2, 

notice that there are only edges i ^ i + 1 and i — > z — 1 in CJ, and that the cycles 
{N + l,N + 1 + 1) and (-iV -l,-N -l-l) belong to Sat($) for any I > 0. 

If A: > then all cycles {i, i + 1) ioi i > N and (i, i — 1) for i < —N have weight 
— fc < and they should not be in Sat($). In particular, Sat($) does not contain 
the edge iV' — > iV' + 1, since it contains N' + 1 ^ N' . However, N' has an outgoing 
edge in Sat($), which must be N' ^ N' - 1. Then N' — 1 ^ N' cannot be there if 
iV' - 1 > A^, so Sat(<I>) contains the edge A^' - 1 ^ A^' - 2 as well. Proceeding this 
way we obtain that Sat(<I>) contains all backward edges i + l^iiorN<i< N'. 
Similarly Sat($) contains all forward edges i — 1 — !■ z for —A" < i < —N. Since A^' 
and A^" are arbitrarily large, part 1. follows. We also obtain that in the restriction 
of G to [— A^, A"] and more generally to [—N — i,N + i] where i > Q there are no edges 
of Sat($) pointing outside of the interval, hence part 2. Indeed, by the outgoing 
edge property, for each node in [— A— i, N+i] there is an outgoing edge, which has to 
point, by part 1., to another node in [— A^ — z, A"+z]. For part 3, Theorem I 1 . 1 1 implies 
that the greatest total weight of a two-cycle does not exceed zero. It amounts to 
show that Sat(<I>) contains cycles, which necessarily have zero total weight. This 
follows from part 2, since the restriction of Sat($) to [— A^ — 1, A + 1] is finite and 
each node has an outgoing edge. Hence the greatest total weight is zero, and (|2.6p 
follows from ((2^ . □ 

Part 1. shows that (A$) implies the bound state condition when fc > 0. More 
precisely, it implies that = ^i+i for all i > A^ and ^i-i = ^i — k for all i < —N. 
When fc > 0, the condition Z > 1 in part 2. can be replaced by Z > 0. In the same 
vein [—A — 1, A^ + 1] can be replaced with [— A^, A^] in the statements below, when 
fc > 0. 

Denote by ^[Ni,N2] the restriction of $ to the interval [A^i, A'2], and by A^j^^ j^^] 
the submatrix extracted from the nodes in the interval [Ai , A'2] . 

Proposition 2.7. Suppose that U satisfies assumption (AU) and fc > 0. Then 

1. If u satisfies A[-7v-i,Ar+i] ^ v = v, then it can be uniquely extended to $ 
satisfying (A<I>) and A ® <I> = $, such that ^[-Ar-i^Tv+i] = v. 

2. If <i> satisfies A®<i> — $, then <i>[_7v-i,Ar+i] satisfies >l[_7v-i,7v+i]®*i'[-jv-i,jv+i] = 

*i'[-JV-l,7V+l]- 

Proof. 1.: By Proposition 12.61 part 1, every eigenvector satisfying (A$), has to 
follow = $i for i > N and — k = for i < —A. Hence it is uniquely 
determined by $[_Ar_i_Ar_|_i]. We also observe that <i>i+i — implies ^i — k < 
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and that — k = implies $^-1 < $i which makes it possible to extend v, 

satisfying A[_jv-i,Af+i] (8> w = u, to $ which satisfies both yl (g) $ = $ and (A$). 

2: If ^ (g) $ = $ then Af^^ jva] ^ ^[Ni.N2] < ^IWi.iVs] for any Ni,N2. By Propo- 
sition part 2, for each i £ [-N - 1, + 1] there is j G [-N -1,N + 1] such 
that Aij + $j = implying that A[_i^_i^,m^i] (g) ^[-n-i,n+i] > ^[-n-i,n+i]- 
Combining with the reverse inequality, we obtain part 2. □ 

For a (possibly infinite-dimensional) matrix A, the Kleene star is introduced as 
in dm]) by 

A* -.^ I ® A® ® . . . . 

In the infinite-dimensional case, this may have infinite number of terms. However, 
in our case the number of terms is always finite for any entry of the Kleene star, 
and the weight of the entry A*j equals to the greatest total weight (i.e., sum of 
weights of the edges) among all paths connecting i to j. 

As in the introduction, by the critical columns of A* we understand the columns 
of A* with indices taken from the critical graph of A. 

Theorem 2.8. Let U satisfy (AU). The set of eigenvectors A ® ^ — ^ satisfying 
(A<^) is nonempty if and only if k is given by (|2.6p . In this case it is the set of 
max-linear combinations of the critical columns of A* , which can be also computed 
as unique extensions of max-linear combinations with the same coefficients, of the 
columns o/ (74[_jv-i,7V+i])* with the same indices. 

Proof. First note that (v4[_7v-i,7V+i])* = (^*)[-Ar+i,w+i] • Indeed, for i,j e [— A^ — 
1,N + 1], if a path connecting i to j has nodes outside [— A^ — 1, A^ -I- 1] then 
it contains cycles. These cycles can be cancelled preserving connectivity of the 
path and not decreasing its total weight, until all nodes of the path are in [— A^ — 
\,N + 1]. Thus, for any such i,j there exists an optimal path which is entirely in 
[-N - 1, A^ -I- 1]. Then (A[_Ar_i,Ar+i])* = {A*)[_]s+i,n+i] follows entrywise by the 
path interpretation of Kleene star. 

Next, it can be verified that any column of A* with an index in the critical graph 
satisfies (A$), and then so does any max-linear combination of these columns. 

Any eigenvector v of ^[-at-i.at+i] is a max-linear combination of the columns 
of (^[-AT-i.Ar+i])* with indices in the critical graph, see Theorem 11.21 Since 
(j4[_jv-i,Ar-f 1])* = (^*)[_Ar+i_Ar+i] , the max-lincar combination with the same co- 
efficients of the columns of A* with the same indices, is an extension of v satisfying 
(A<I>). Such extension is unique by Proposition 12.71 part 1. Conversely, by Propo- 
sition 12.71 part 2, any vector satisfying v4 O $ = $ is a unique extension of an 
eigenvector of 74[_jv-i.Af+i] • D 

2.3. Undressing transform. In this subsection, assumptions (A$) and (AJ7) are 
assumed everywhere. 

It is easy to see that each soliton corresponds to a critical component in the 
associated digraphs of ^(7) and A(6). Hence it follows that each soliton gives rise 
to a pair of fundamental eigenvectors of ^(7) and A{5) which we denote by <I>(^) 
and We next examine the transformation U 1-^ U defined by 

(2.8) ^^:^^^ + ci,W +$f)-ci>«-<)^. 

To compute this transformation explicitly, we need the relations between neigh- 
bouring coordinates of ^^^^ and They are as follows. 
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Proposition 2.9. Suppose that (CI) holds. Let (/,..., / + s) be a sohton and let 
and <i>(^' be the pair of fundamental eigenvectors of ^(7), resp. A{6) associated 
with it. Then 



*f ^ - *!+i = -k, t<l + s-l, 
*!+i - =u„ t>l + s-l. 



Proposition 2.10. Suppose that (C2) holds. Let {I, . . . ,1 + s) be a soliton and let 
<i>(^' and <i>*^^) be the pair of fundamental eigenvectors of ^(7), resp. A{S) associated 
with it. Then 



Ui - 1, if Ui + Ui^i < 1, 



(2.10) 



if 


Ui - 




-1 > 1 


if 


Ui - 




-1 < 1: 


if 






-1 > 1 



"^+1' ' -^-^ *>z, 



f2i ci) u,_i — 1, if Uj + < 1, 

y-Ui, it + Ui_i > 1 



1 - Ui+1, if Ui + Ui+i > 1 



Proof. (Propositions 12.91 and I2.10p In both cases, we essentially have to examine 
which edges of ^(7) and A{5) are in the saturation graphs of and Then 
we use the explicit formulas for the coefficients of ^(7) and A{5), see Lemmas 12.11 
andO 

In both cases (CI) and (C2), the cycle {I, l + l) is critical in ^(7), thus can 
be chosen as Z or Z + 1 column of ^(7)*. If J < ' then 



(2.11) -^f - <i>Ji\ = A(7)*,+i = ^(7)m+i- 
Hence (^''p - "I>|^\ = ^(7)^.^+1 for all i < I. Analogously for j > I 

(2.12) <i>f ^ - = A(7)*, = ^ A(7)m-i. 

i=l+l 

Hence - = ^(7)^+l,^ for all i > I. 

In both cases (CI) and (C2), the cycle {I + s — 1,1 + s) is critical in A{5), thus 
^'^^^ can be chosen as Z + s — 1 or Z + s column of A{S)*. Arguing as above, we obtain 
that $f' - = ^('5)»,»+i for alH < / + s - 1, and that - $f' = ^(^)^+l,^ 
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$(1) - 


^i+l 


= A(7),,i+i, 


i < I, 


^«+l 






i > I, 


$f ) - 


-$!?i 




i <l + s 


'j>i?i 


2 




i>l + s 



for all j > / + s — 1. Summarizing we have: 



(2.13) 



It remains to use the exphcit expressions for coefficients of ^(7) and A{S). □ 

Next we estabhsh exphcit expressions for undressing (|2.8p . in the situations (CI) 
and (C2). 

Theorem 2.11. Let (Z, ...,/ + s) be a soliton and let and <i>(2' be the pair of 
fundamental eigenvectors 0/^(7), resp. A{5) associated with it. Then in the case 
of (CI) 

^fi<l 

(2.14) Ui = ( Ui^i — Ui+i, if I < i < I + s, 




if i > I + s, 



and in the case of (C2) 



Ut-i, ifi<l, 
(2.15) Ui=<(l-Mj, ifl<i<l + s, 

^Ui+i, ifi>l + s. 

Proof. The computation is straightforward, using ()2.8p and (j2.9p in the case of 
(CI), or (1211) and (l2J0ll in the case of (C2). 

Namely in the case i < I we use the first and the third relations of (|2.9p 
and (|2.10l) . substituting them into (|2.8p . In the case i > I + s we use the sec- 
ond and the fourth relations of p.9p and (|2.10l) . 

In the case I < i < I + s, we use the second and the third relations of ()2.9p 
and (|2.10p . If (CI) holds, note that for / < i < Z -f- s the computation yields 
Ui^i + Ui-i + Ui ^ k. However all cycles {i,i — 1) and {i,i + 1) are critical with 
+ Ui — k = Ui + Ui-i — fc = 0, hence we obtain = Wi+i as in (|2.14p . If 
(C2) holds, then the computation (use the case Ui + Wi+i > 1 since we are inside 
the soliton) yields 2 — li^ — fc, which is 1 — since k — 1. □ 

We obtain that in the case of (CI), according to (|2.14p . the selected soliton 
{I, . . . ,l + s) loses two units of its length, and the rest of the potential U gets shifted 
to the right before the selected soliton, and to the left after the selected soliton. 
The remaining part of the selected soliton also gets shifted, but the direction does 
not matter since = Ui+i for aW I < i < I + s. 

In the case of (C2), according to ()2.15p . the selected soliton also loses at least 
2 units of length on the ends, and the remaining part may shrink and brake into 
several solitons. Indeed, we have 

ui + iii+i = 1 - ui+i + ui-i <ui+ ui-i < 1 

and also 

Ul+s-l + Ui+s = 1 - Ui+s-l + Ul+s+1 < Ui+s + Ui+s+1 < 1- 
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Ul+l Ui+2 Ui+3 Ui+4 

ui-i — k ui — k it;+i — k — k 




Sat($(2)) : 1 



Ul+l Ul+2 Ul+Z 

U: ui Ul+l = ui+2 = ui ui+3 ui+i 



Figure 4. Undressing in (CI) 



Sat($(i)) 



Sat($(2)) 



U: 
U: 




1 - Mi 1 - Ul+l 1 - Ul+2 

Ul^l - 1 -Ul + l -Ul+2 



1 - ";+2 1 - Ui+3 

Ul Ul + l Ul+2 Ul+3 

Ul-l 1 - tli+i 1 - Ul+2 ""(+4 



Ul+i 
-Ul+3 



Ul+3 



Ul+4 
Ul+5 



Figure 5. Undressing in (C2). The dashed edges indicate that 
they belong to saturation graphs if and only if Ui + Ui+i = 1 for 
the corresponding i, that is, if the two-cycle containing them is 
critical. 



For i in I < i < I + s — I we obtain 

Ul + Ul+l = 2-Ui- Ul+l, 

which is not less than 1 only if Ui + ui+i — 1. 

An important special case of solitons in case (C2) is when they are of the form 
(a 1 ... 1 b) with a, 6 > 0, and when all elements between them and outside the 
soliton area equal 0. In this case, such soliton completely disappears turning into 
background after the corresponding undressing transform (which justifies the name 
"undressing"). 

We give a graphical example of undressing where the selected soliton has length 
3. Figure m demonstrates undressing in the case (CI): look at the difference between 
U and U . Figure [5] demonstrates undressing in the case (C2). 



AN APPLICATION OF THE MAX-PLUS SPECTRAL THEORY 



13 



3. Adding constraints 

In this section we verify whether a fundamental pair also satisfies the 

last two equations of (|l.ip . which we rewrite as 

(3.1) <i>W = max(<i>[^\, <i>W +ui- 1), 

(3.2) = max($[^^ - a, $f ^ + ui + k - 1). 

The parameter cr, equal to — fc in (jl.ip . will be specified later. Let us remark so 
far, that since = = at all large enough I, using p.2p we obtain cr > 0. 

3.1. Case (CI). 

Proposition 3.1. In case (CI), any fundamental eigenpair satisfies p.ip and 



Proof. In the case (CI), the graph Q{S) is the same as the graph Q{-f) shifted one 
position to the right, and the same is true about the graphs Sat($'^') and Sat($'^') 
for any fundamental eigenpair $(2). Then — = $|^2 ~ ^i+i for all i 
and hence — for all i. We put cr = and verify the remaining inequalities 



(3.3) 



$!+\ > $f ' +u, + k-l 



comparing them with p.9p . The verification follows from k < 1 and Ui + u^+i < 
1. □ 

Thus in the case (CI) any fundamental eigenpair satisfies (jl.l[) . Also note that 
due to max-plus linearity, any max-plus combination of fundamental eigenpairs 
($(1)^ $(2)) ig again a solution of the system, so that in general the solution space is 
highly degenerated. Note that this result justifies the study of undressing by means 
of fundamental pairs in the case (CI). 

3.2. Case (C2), one soliton. Here we verify that in the case when there is just 
one soliton in case (C2), (/,..., Z + s), the fundamental pair satisfies (|3.ip and (|3.2p . 

In the following table, we consider an example where the potential (i.e., so- 
lution of udKdV) consists of one soliton (tti 7r2 tts 714), where the real numbers 
TTl , 7r2 , 7r3 , 7r4 < 1 are real numbers such that 7ri-|-7r2 > 1, 7r2+7r3 > 1 and 7:3+774 > 1. 





I : 







1 


2 


3 


4 


5 


6 




Ul : 







TTl 


7r2 




7r4 








11 


— min(u;, 1 — 


ui-i) : 





TTl 


1-^1 


1 - 7r2 


1 - TTs 








Si 


= min(u;_i, 1 


- Ul) ■■ 








1 - 7r2 


1 - TTs 


1 — 7r4 


71'4 






The digraphs ^(7) and G{6) are displayed on Figured 

The saturation graph Sat(<l>(^^) is a subgraph of ^^(7) shown in the upper part 
of the picture: the edges not belonging to it are drawn as dashed. Similarly, the 
saturation graph Sat($(2)^ jg g, subgraph of G{S) shown in the upper part of the 
picture. Combining these graphs, we see that we have to analyse three cases: 1) 
before the soliton area (to the left), 2) in the soliton area, 3) after the soliton area 
(to the right), with three different combinations of edges used by the fundamental 
pair. 
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TTl 1^- TTl 1 - 7r2 1 - TTs \ 

before soliton \ soliton area \ after soliton 

^1 ^1 \ —7^2 ^T^S —n^'- 77^ — 1 




0: 



71"! 



7r2 



7^3 



774 



Figure 6. The case of one soliton 



To formalize the arguments let us introduce the notion of soliton area in general. 
Let (Z, ...,/ + s) be a soliton. By the soliton area we mean a pair of subgraphs of 
G{'j) and Q{S): 1) the subgraph of G{j) extracted from the nodes I, . . . ,1 + s, 2) 
the subgraph of Q{S) extracted from the nodes I + 1, . . . ,1 + s + 1. 

Theorem 3.2. // U contains just one soliton (I, ... ,1 + s), then the pair of funda- 
mental eigenvectors associated with it satisfies (|3.ip and (|3.2|) . 



Proof. We are going to use relations between the neighbouring components of 
and written above in (j2.10|) . 

For the area after the soliton we obtain 



(3.4) 



(1) 



and in particular ^ -^i-^f ■* = *J'l+2~*J'i+i fo^' * ^ This implies $f ■* = 

follows from 'l',-^^ > ^'P' — ct, and 



We also obtain 



for the area after the soliton. Equation 

(2) (2) 

since = $j + Ui by the second equation of 
(3.5) <i>P' - $W > 1 

from the first equation of (|3.4p . since Mi + u^+i < 1 for i > I + s, which together 
with $f ^ = makes (jXT]) . 
For the soliton area we obtain 

(3.6) 



1 - Ui 



and in particular $f - <i>-J^\ > <i>-^\ - $-^2. Equation ([3?6|) implies that 



.(2) 



(2) 



(3.7) 
Defining 

(3.8) 



(2) _ 



*r'-*l+i + (^^«+«'^-i-i)- 



Z+s-l 

cr = ^ (u; + - 1) 
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we obtain that 

(3.9) $p'-a<a>!+\<<i'r^ 

To show p.2p and (13. ip we observe that (13. 6p furnish the remaining necessary 
equalities = $f ^ + u, and $f ^ = + - 1. 
Before the soliton area we obtain 

(3.10) 



$(2) _$(2) ^ ^ 

Here the difference '^^p — is stable, and by p.7p and p.8p it equals to cr, so 

= $1^' - cr for all i < I. Equation I^^Q follows from $f ^ > and since 

(1) 

i+l 



^f'^ = + - 1 by the first equation of (I3.10p . We also obtain 
(3.11) $!+\>$r+ti« 

from the second equation of p. 101) . since Ui + Ui-i < 1 for i < I, which together 
with = - cr makes □ 

This result implies that in the case of one massive soliton, when U satisfies (At/), 
a solution $ to (|l.ip satisfying (A$) exists and is unique. 

3.3. Case (C2), several solitons. We have seen above that in the case of one 
soliton, the last two equations of (II. ip are satisfied automatically. However, the 
graphs ^(7) and Q{5) contain edges which are dangerous to use. If Sat($(^)) or 
Sat($(^)) contain such edges then the last two equations of ()l.ip are violated. 

Lemma 3.3. Let ^l^^^ and be a solution to (ITTT|) . Then Sat($(i)) cannot 
contain edges (i, i + 1) if Ui + Ui^i > 1, and Sat($('^^) cannot contain edges (i + 1, i) 
if Ui + Ui+i > 1. 



Proof. We use Lemma [2.21 being in the case of (C2). 

If Sat(3>(^)) uses + whenwj + u,_i > 1 then $f ^ - = By (f3T|) 

we should have ^[^"^ — > Ui — 1 and hence Ui + < 1, a contradiction. 

IfSat($(2)) uses (i + l, i) when Ui+Ui+i > 1 then -$f ^ = l-u^+l. By 

(2) (2) 

we should have — ^1 > ui and hence Ui + Ui+i < 1, a contradiction. □ 

It can be checked that the use of other edges does not lead to such contradictions, 
and also, using Lemma [^Tl that there are no contradictions in the case of (CI). 
The following negative result is now easy to see. 

Theorem 3.4. Let U satisfy (C2) and contain more than one soliton. Then no 
pair of fundamental eigenvectors associated with a soliton can be a solution of (|l.ip . 

Proof. According to Lemma 13. 3[ to each soliton there corresponds a number of 
consecutive forward edges in ^^(7) that cannot be used by Sat($^^^), located im- 
mediately after the corresponding critical cycle in ^(7). Further, there is also a 
number of consecutive backward edges in Q{S) that cannot be used by Sat(<i>(^)), 
located before the corresponding critical cycle in G{S). 
If $(1) and $(2) 

are a pair of fundamental eigenvectors, then the switch from 
backward to forward edges can happen only once. To avoid all forbidden forward 
edges after the critical cycle of ^(7) corresponding to the first (i.e., left- most) 
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— 1 TTl — 1 — TTl —1 —1 TTs — 1 — TTa —1 




TTl 1 - TTl 

soli 



TTa 1 - TTs 

sol2 




0: 



71"2 



7r4 



Figure 7. The case of two solitons: "dangerous" edges 

soliton, Sat($'^^^) has to use backward edges only, which implies that (l>(i),$(2) 
should be the pair associated with the first soliton. But then also Sat(<i>^^-') uses all 
backward edges after the first critical cycle oiQ{5), including all forbidden backward 
edges corresponding to the next solitons. □ 

In the following table, we consider an example where the potential (i.e., solution 
of udKdV) consists of two solitons (tti 1:2) and (tts tt^), and the real numbers 
""i, 7i'2, TTs, 7r4 < 1 are such that tti + tt2 > I and 773+774 > 1. 



I: 01 2 345 6 78 

ui : 7ri 773 774 

7; = min(it;, 1 — : 7ri 1 — 7ri 773 1 — tt^ 

(5; = min(ui_i, 1 — u/) : 1 ~ tt2 '^2 1 — 774 774 

Digraphs Q{-f) and Q{S) are displayed on Figure [71 

Here the bold edges are used by the pair of fundamental eigenvectors corre- 
sponding to the first soliton {tti, 712)- The red edges are the ones which lead to a 
contradiction with (|3.ip (graph ^^(7), upper part of the figure) or (|3.2p (graph Q{5), 
lower part of the figure) . We see that the pair of fundamental eigenvectors has to 
use one of the forbidden edges, hence it cannot satisfy both (13.11) and p.2p . 

4. Conclusions and projects 

In this paper we attempted to build the max-plus theory of (jl.ip . Based on the 
observation that the first two equations represent max-plus spectral problems, we 
explained how the finite-dimensional max-plus spectral theory applies to them. We 
studied pairs of fundamental eigenvectors associated with each soliton, describing 
the undressing transform and showing that these pairs yield a solution of (jl.ll) in 
some situations. 

The remaining nontrivial case is when J/*^*-* has several massive solitons, where 
we have shown that the pairs of fundamental eigenvectors violate the last two 
equations of (jl.ll) . Willox et al. [TT] report that a solution can be found also in this 
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nontrivial case. It is desirable to work out a systematic comprehensive approach to 
solving in this case, and in particular, to understand whether a max-plus linear 
combination of fundamental pairs could be a solution. Then one could proceed with 
the study of undressing transform associated with any solution of (|l.ip . and the 
details of application of this theory to solving the ultradiscrete KdV equation (|1.2I) . 
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